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Abstract The two-stream instability is a common instability, responsible for many observed phe- 
nomena in nature, especially the interaction of jets of various origin with the background plasma 
(e.g. extragalactic jet interacting with the cosmic background). The dispersion relation not consid- 
ering magnetic field is described by the well-known Buneman relation. In 2011, Bohata, Bfen, and 
Kulhanek derived the relation for the two-stream instability without the cold limit, with the general 
orientation of magnetic field, and arbitrary stream directions. The maximum value of the imaginary 
part of the individual dispersion branches uJ n (k) is of interest from the physical point of view. It 
represents the instability growth rate which is responsible for the turbulence mode onset and sub- 
sequent reconnection on the ion radius scale accompanied by a strong plasma thermalization. The 
paper presented here is focused on the non-relativistic instability growth rate and its dependence on 
various input parameters, such as magnitude and direction of magnetic field, sound velocity, plasma 
frequency of the jet and direction of the wave vector during the jet - background interaction. The 
results are presented in well-arranged plots and can be used for determination of the plasma pa- 
rameter values close to which the strong energy transfer and thermalization between the jet and the 
background plasma occur. 
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The most common plasma instabilities are the two-stream instabilities, which can occur during a plasma jet 
interaction with the plasma background. Such situations are observed in astrophysical processes, e.g. interaction 
of the galactic jets with the intergalactic medium (e.g. Silketal. 2012 and references therein) or interaction of 
star jets with the interstellar medium ( Murp hy" et al. 2008| l. Oskar Buneman derived the basic dispersion relation 
describing such instabilities in the late 1950's for cold unmagnetized plasmas (Bunem an 1959b . The magnetohy- 
drodynamic instabilities in the ideal plasma are discussed in ( |Bonanno & Urpin 201 1) . Magnetic fields are crucial 
for the phenomena taking place in jets ( |Urpin 2006] l. In 2011, Bohata et al. published a paper containing the 
derivation of the non-relativistic dispersion relation for magnetized plasmas and without the cold limit restriction 
(Boha ta"et al. 201 lb . It is called it the Generalized Buneman Dispersion Relation (GBDR) and it is described by 
the equation 
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is the Doppler shifted frequency, ui ca is the cyclotron frequency, ui pa is the plasma 

frequency, Fa' is the Lorentz force, eB is the unit vector in the direction of the magnetic field and c sa is the sound 
velocity. Index a denotes the corresponding media (jet or background). 

In the previous work (Bohata et al. 201 1), the numerical solution for the case of two identical plasma beams 
with the same velocities, but opposite directions was found for various input parameters. The situation of the 
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plasma jet penetration into the plasma background was studied as well and the numerical solution for this problem 
was found Porky 20 1 2a} 120 1 2bb . 

The maximum of the imaginary part of the solution is denoted as the Plasma Instability Growth Rate (PIGR). 
This paper is focused on finding the plasma parameters for which this maximum occurs (the non-relativistic case 
and the plasma jet interaction with the plasma background are assumed). The calculations are performed on the 
microscopic level and using the linear approximation. For the plasma parameters leading to the maximum of the 
imaginary part of the dispersion relation, the instability arises and amplitudes of all variables grow exponentially. In 
such a situation, the linear approximation is no longer valid, and other methods for modelling of the physical phe- 
nomena must be introduced. One of the possibilities is the Particle in Cell (PIC) simulation, (e.g. Stockem 2008). 
The results of these calculations can be therefore applied to: 1) The search for the instability regimes in which 
strong thermalization, turbulence, micro-reconnections on ion radius, non-thermal radiation, shock onset and other 
interesting phenomena can occur. The subsequent PIC simulations of the plasma behaviour leading to significant 
phenomena seems to be the most reasonable next step of the research in this regime. 2) The tests of the accep- 
tance of the PIC codes (the PIC code must lead to an instability onset for the parameters calculated by the method 
proposed in next paragraph). 

1 METHOD 



Rhe following indices were designated in our analysis: "j" for parameters of the jet and "b" for parameters of the 
background. It is beneficial to transform the variables and the whole GBDR relation to the dimensionless form. 
After this step the relations simply scalable and the equations are covariant against to this transformation. It suggest 
that the results can be used for both space and laboratory plasmas, such as thermalization in astrophysical jets or 
in fusion experiments. The relations for the dimensionless form were chosen with regard to the zero background 
velocity as ( |Horky 2012all2012bb : 
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J?j = uj — k cos tp sin 0k , 
J?b = TD — ku 2 cos tp sin # k . 

The reference system was set according to the Fig. [I] in which the directions of the respective vectors u Q , B and k 
are drawn. The wave-vector can point in any direction, the magnetic field vector lies in the (x-z) plane and the jet 
is directed along the jt-axis. The vector coordinates are 

u Q = (u a ,0, 0) , 

B = (5 sin B) O,.B cos B ), (3) 
k = (k cos if sin 0k, k sin tp sin 0k, k cos 0k) ■ 
After simple manipulation the dimensionless form of the dispersion relation becomes (Horky 2012a, 2 1 2bl > 
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where the goniometrical terms were denoted as 

G\ = (cos 0b sin ip sin 0k) , 
G 2 — (cos ip sin k sin B + cos 0k cos 0b ) , 
G3 = (cos 2 0b cos ip sin 0k — cos 0b cos 0k sin 0b) 



(5) 
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Fig. 1 The system of coordinates used in calculations. 

This dimensionless relation is a polynomial equation with complex roots of the 8 th order. The algorithm developed 
by Hubbard, Shleicher, and Sutherland (Hubba rd et al. 200Tb was used to find the solution. The algorithm was 
implemented in the Wolfram Mathematica 8.0.1. and in its principle ensures that the iteration from the initial seeds 
will converge tol the roots in comparison with much more simpler Newton-Raphson method. The results were 
arranged into plots in which the real branches of the solution have different colour than the imaginary branches, 
and the maximum imaginary value determining the PIGR value was highlighted. The example of the program 
output is shown in the Fig. [2] 

In the next step, the PIGR value dependence on various parameters of the dimensionless GBDR (such as 
cyclotron frequencies of the jet and the background, the sound velocities of the jet and the background, the plasma 
frequency of the jet and the directions of the magnetic field and of the wave-vector) was found. 

2 RESULTS 

The PIGR value was calculated during the program cycle running from the minimum to the maximum value of the 
tracked parameter while other parameters were fixed at their initial values. Intervals of these parameters are shown 
in Tab.[T] It was not necessary to change the jet velocity, because its dimensionless value was fixed at 1. 



Parameter 


Initial value 


Minimum value 


Maximum value 




0.5 


0.1 


3.0 




0.1 


0.1 


1.5 


^Pj 


1 


1 


5 


0k 


tt/2 





tt/2 


¥> 








tt/2 


6b 


tt/4 





tt/2 



Table 1 Parameters used for the numerical solution. The dispersion relation for the initial values is 
depicted on the Fig. 1. 



2.1 The PIGR value dependence on the cyclotron frequencies 

The cyclotron frequency of the jet and the cyclotron frequency of the background were increased from the min- 
imum value of 0.1 to the final value of 3.0 with the step size 0.1. The cyclotron frequency is proportional to 
the magnetic field intensity influencing the charged particles. The PIGR value dependence on the cyclotron fre- 
quency of the jet is depicted in the Fig. [3] where the almost almost linearly increasing character of this relation for 
uJ C j > 0.6 is noticeable. The change of the slope at this point (ZJ C j = 0.6) corresponds to the localization of the 
minimum of the two different imaginary branches of the dispersion relation. The PIGR value dependence on the 
cyclotron frequency of the background is more complicated than the case of jet cyclotron frequency. In the Fig. |4] 
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Fig. 2 Real (blue) and imaginary (red) branches of the GBDR dispersion relation and PIGR value 
(denoted as Max) for uJ C j = <U c b = 0.5, c S j = c s b = 0.1, w p j =1,8^ = ir/2, ip = 0, Ob = 7r/4. These 
values were used as initial values for the calculations, see Table 1 for details. 



the descending character for lower frequency values is visible. The curve reaches minimum and then it rises to an 
asymptote. The minimum was numerically determined to be uJ C 6 = 1.313 and the corresponding PIGR value was 
equal to 0.39795. This effect is caused by the fact that the solution has two imaginary branches in this area and 
while the value of uj c i, is increasing the first branch is descending and the second is rising. In the minimum both 
branches have equal PIGR values. 

2.2 The PIGR value dependence on the sound velocities 

Sound velocity is proportional to (Tq/toq,) 1 / 2 , where T a is the plasma temperature, and m a is the mass of the jet 
or of the background particles (electrons or ions). The index a labels the corresponding media (jet or background). 
Modification of the original Buneman dispersion relation by addition of the sound velocities of both media is 
a result of the calculation with non-zero pressure, i.e. without the cold limit. The dimensionless parameter c s 
involves the plasma jet velocity, see Eq. (|2jl, and c S j > 1 indicates subsonic jet and vice versa. Both the sound 
velocity of the jet and the sound velocity of the background were increased from the initial value 0.1 to the final 
value 1.5 with the step 0.1. The PIGR value dependence on the sound velocities of both jet and background is 
depicted in the Fig. [5] The jet dependence (red circles) has descending character and PIGR value is zero, while 
c S j > 1. It implies that for a subsonic jet (in dimensionless form the sound velocity equals 1) the GBDR relation has 
no imaginary branch and therefore the PIGR value is zero and no instabilities occur. The PIGR value dependence 
on the sound velocity of the background (blue circles) is more complicated. The interesting peak resides at the 
value c s b = 0.9. We made three-dimensional plot of the imaginary branches of the GBDR solution to uncover 
the origin of this local maximum. The result can be seen in the Fig. [6] First axis corresponds to c s b, second to k, 
and the vertical axis to the value of the imaginary branch of the PIGR coefficient. It is clear now that the peak 
originates from the ridge present at the dispersion relation solution. 
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Fig. 3 The PIGR value dependence on the jet cyclotron frequency 
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Fig. 4 The PIGR value dependence on the background cyclotron frequency 
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2.3 The PIGR value dependence on the plasma frequency of the jet 

All dimensionless frequencies in the system are related to the background plasma frequency, see Eq. Q. It means 
that the dimensionless plasma frequency of the background tZ7 p b is equal by definition 1 and dimensionless plasma 
frequency of the jet cj p j is in fact the ratio of the jet and the background plasma frequencies. This parameter is 
therefore proportional to the (riej/rieb) 1 ^ 2 - During numerical calculation it was increased from the initial value 1 
to the final value 5 with the step size 0.5. It is a rather big step, but as it can be seen i then Fig. [7] the dependence 
is very simple and without any discontinuities or local maxima or minima. 



2.4 The directional PIGR value dependencies 

The PIGR value dependence on the magnetic field direction is simply predictable from the Lorentz equation of 
motion. Longitudinal magnetic field will evoke less disturbances than the perpendicular one. As you can see in the 
Fig. [8] the PIGR value has maximum at 6b — (perpendicular direction) and decreases for increasing 6>b- The 
PIGR value dependence on the direction of the wave vector is also predictable due to the dot product between k 
and u Q in the GBDR relation, so if the angle between the wave vector and the velocity equals 90°, the PIGR value 
should be zero. In the Fig. [8] the dependence has descending character and it's zero at the angle 90°. Because of 
the cylindrical symmetry, both angles (p^ and 6>b were altered only from 0° to 90° with the step of 10°. 
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Fig. 6 The imaginary branches dependence on c s b 
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Fig. 7 The PIGR value dependence on uJ p j 
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Fig. 8 The PIGR dependence on the magnetic field direction (blue circles) and on the wave vector 
direction (red circles) 



3 CONCLUSION 



Plasma jets from black holes and other types of astronomical objects are driven by the magnetic fields, and 
classical Buneman instability analysis (without magnetic fields) is inapplicable. All calculations must be per- 
formed using the Generalized Buneman Dispersion Relation (GBDR) with nonzero pressure and nonzero mag- 
netic field. The Plasma Instability Growth Rate (PIGR) as the maximum of the imaginary parts of the GBDR 
relation was numerically calculated in this paper. The PIGR value is responsible for strong thermalization during 
the jet-background interaction and these calculations can be useful for understanding of the underlying processes. 
Furthermore, the known PIGR value can be used as a simple test of Particle in Cell (PIC) numerical methods fre- 
quently used for the plasma jet simulations. It is an interesting but still open question if the PIGR value could be 
calculated analytically directly from the dimensionless GBDR relation. The dispersion relation is not anisotropic 
in the velocity space. This possibility can cause other phenomena, e.g. particle acceleration, shock origin, etc. 
(Nishikaw a et al. 200 6; Mizun oet al. 2009 and references thereinl l. which will be the topic of detailed PIC simu- 
lations. 
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